Topology, mass and Casimir energy 

N. Ahmadi"* and M. Nouri-Zonoz ^ 

Department of Physics, University of Tehran, 
North Karegar Ave., Tehran 14395-547, Iran 
^ Institute for studies in theoretical physics and mathematics, 
P Box 19395-5531 Tehran, Iran. 

Abstract 

The vacuum expectation value of the stress energy tensor for a massive scalar field with arbitrary 
coupling in flat spaces with non-trivial topology is discussed. We calculate the Casimir energy in 
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evaluation of the Casimir energy consists in an expansion in terms of the lengths of these paths. 
We will show how different paths with corresponding weight factors contribute in the calculation. 
The optical approach is also used to find the mass and temperature dependence of the Casimir 
energy in a cavity and it is shown that the massive fields cannot be neglected in high and low 
temperature regimes. The same approach is applied to twisted as well as spinor fields and the 
results are compared with those in the literature. 
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I. INTRODUCTION 



Investigations of the different aspects of quantum field theory in spaces with non-trivial 
topology I] may be interesting for several reasons. First, it is a possibility that the early 
Universe was in fact multiply connected, so the combined effect of topology and curvature 
is likely to have been of considerable importance in the early stages of the evolution of 
the universe. Second, in theories with extra dimensions such as the well known Kaluza- 
Klein theory, extra dimensions may need to be compactified providing the space with non- 
trivial topology. Furthermore field theory in spaces with non-trivial topology (the one- 
torus compactified space) is quite similar to the field theory at non-zero temperature. The 
Casimir effect on spaces with non-trivial topology has been investigated in different contexts 
(Q] and references therein). Most of the previous studies are concerned with one compact 
dimension (representing imaginary time or a compact spatial dimension) in either a curved 
or flat spacetime Oj. The issue of obtaining the right order of magnitude for the possible 
contribution of the Casimir energy density to the cosmological constant, corresponding to a 
free scalar field, is considered in toroidal universes with p large and q small dimensions jj] . 
Since the Vacuum polarization appears in manifolds with non-trivial topologies as well as 
in bounded domains with boundary conditions, the question of the definition of the vacuum 
energy is a most fundamental one in any quantum field theory. On the other hand imposing 
(mathematical) boundary conditions on physical quantum fields turns out to be a highly 
non-trivial act. The difference between the zero-point energies corresponding to the scalar 
field when boundary condition is enforced and that in the absence of any boundary condition, 
diverges as the boundary is approached. Integrating over all space yields an infinite vacuum 
energy per unit area of the boundary surface. Since the infinite terms do not depend on the 
distance between the plates, these divergences do not usually show up in the calculation of 
Casimir force. To obtain the finite results through the use of a suitable regularization, many 
efforts have been done to identify and separate the singularities by physically meaningful 
cutoffs j^. The whole argument is based on boundaries that constrain all modes of the 
fluctuating field, introducing divergences into physical observables. Therefore one needs a 
realistic boundary condition through a high energy cutoff that limits the fluctuations and 
thereby introduces a high energy scale into the theory. Modes with energy below this scale 
obey a boundary condition while modes with energies at or above this scale do not. Hence 
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the boundary condition applies to the physics at energies much lower than this scale where 
there is no need to reference the scale, as it was the case in the Casimir's original problem. 
In this paper we examine the Casimir-type effects in four dimensional spacetimes with 

n 

R{time) x S (space) topology, using an approach introduced recently 'p\ for an approximate 
calculation of Casimir energy in generic geometries like the Casimir pendulum and a sphere 
opposite a plane. Although an approximation, it leads to exact results for flat geometries. 
It is specially useful when the normal modes cannot be determined. The basic idea is to 
employ the regularization method of point splitting through an expansion of the Casimir 
energy in terms of the classified, closed geodesic paths and then identify those classes that 
contribute to the Casimir effect. The expansion is given in terms of the lengths of the paths 
and the paths with zero-length (corresponding to the high frequency components of the 
field) introduce divergences as they only probe the geometry in the immediate vicinity of 
the spacetime point of interest and in this restricted neighbourhood the topology does not 
change. 

Including a massive field in the Casimir energy calculations is the simplest generalization of 
the usual, massless scalar field calculations. Casimir-type effects in the presence of massive 
fields have been studied in |2] and jS] for both scalar and spinor fields. There it has been 
found that a non-zero field mass always decreases the magnitude of the energy density. Here 
we will encounter cases in which massive fields act in the opposite direction. The outline of 
the paper is as follows, in Sec. II the stress-energy tensor is written in terms of the trace 
of the Hadamard function which is the basic object from which (0|T^j,y|0) is constructed 
by differentiation. We will review the Casimir energy calculations between parallel plates 
as an example of the use of the optical approach in section III. Some applications of this 
approach in fiat four dimensional manifolds, either orientable or non-orientable, are studied 
in section IV. Divergences due to the massive fields at low and high temperature limit are 
also discussed in this section. Studying spinor fields in section V, we show that the optical 
method can yield exact results. Conclusions and summary are discussed in the last section. 
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II. STRESS-ENERGY TENSOR 



The action functional for a scalar field $ in a curved spacetime is 

S [$] = i J ($.^$'^ - ^7^$2 - m'$') d^x (1) 

where g is the determinant of the background metric g^^, TZ is the curvature scalar and <; 
is the gravitational coupling. Varying $ infinitesimally, we obtain the scalar field equations 
[□j; + + <,TZ] $ = 0. The classical scalar stress tensor density is given by 



T, 



;i-3c^)mV^' 



(2) 



2"^ 2 

The transition from classical to quantum fields is made by replacing the classical field $ by 
a field operator. We then note that eq. ^ is constructed from the products of the field 
operators or their derivatives at the same spacetime point. These quantities diverge when 
their vacuum expectation value is taken. To avoid these divergent quantities, one operator in 
each product is moved to a nearby point. The point separated vacuum expectation value of 
the energy-momentum tensor is then expressed in terms of the so-called Hadamard function 
G*^"*"^ (x, x) as, 



(0|T^,|0) = Jlim 

Z x—>x 

1 / 1 3 



3<;) ni^g^^, 



G^^\x,x) (3) 



(4) 



where G*^^^ { X ^ X ^ IS defined by G^^^ {x,x) = (0 |{$ (x) , $ (x)}| 0), and satisfies 

[□^ + m^ + ^7^] G(^) (x, x) = 0. 

Because of the difficulties in calculations of quantum field theory in a general curved space- 
time, it is very natural to deal with some specific spaces. To study the effect of topology on 
Casimir energy, we relax the curvature effect by restricting ourselves to flat spacetimes with 
i? X S topology whose metric is given by 

{dsf = (rfx°)' - {dx'Y - {dx'f - [dx^f. 

For the Minkowski space G« (x,x) is given by Q 

G(^) (x,J) 



(5) 



—Qi2a)KAmy/2^) + ^= 



Q{-2a)h{m^^). 



(6) 



In which O is the step function and Ki and Ii are the modified Bessel functions of the first 
and second kind, respectively. The quantity a (x, x) = \gai3 {x" — x") (x^ — x^) is a biscalar 
equal to one half the square of the geodesic distance between x and x [ill. As discussed in 
the introduction, divergences occurring in the stress-energy tensor calculations are due to 
the high frequency components of the field. In Minkowski space quantum field theory, these 
divergences are simply discarded using the normal ordering prescription for field operators. 
On the other hand one may employ an ultraviolet regulator function e~^l^^ to cut off the 
ultraviolet divergences, and then take the difference between the values of (0 |T^,^| 0) in the 
presence of the boundary or the topological condition and in the absence of them, letting 
A ^ oo at the end of the calculation. Alternatively, in the Green function approach, one first 
subtracts the unbounded Minkowski space expression for G'-^-' (x, x) from the same function 
evaluated in the topology of interest, and then takes the limit x — > x. In this paper using 
the latter approach we expand Green functions in terms of all geodesic paths from x to x, 
classified and labeled by the family index n through their endpoints x„. The lengths of these 
paths contribute with equal or different weight factors, M„, 

G«(x,x) = J]m„G«(x,x„)- (7) 

n 

In this approach, it will become clear that the divergences arise due to arbitrarily short 
length paths which could be isolated and discarded from the outset, leading to the infinity- 
free Casimir energy calculations. 

III. CASIMIR ENERGY BETWEEN THE PARALLEL PLATES 
A. Ideal parallel conductors 

In this section, using the optical approach, we review the Casimir energy corresponding 
to a massive scalar field in the most simple configuration i.e. a rectangular parallelepiped 
with ideal walls of dimension a, h and c with a -C 6, c. We employ Dirichlet condition at the 
boundaries i.e., 

$ (x°, 0, x^ x^) = $ (x°, a, x^ x^) = (8) 
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The stress-energy tensor is given by 



T, 



(1-2^)$;^$;, 



2^ 



9 ,ot^ 



(9) 



Using the field equations along with eq. (jH)), vacuum expectation value of T^^, in a flat 
spacetime could be expanded in terms of geodesic paths from x to 



(0|T^,|0) = -lim 



[1 - 2^) V^V. + ( 2c; - ^ ) g^^VaV' 



(10) 



Every geodesic path contributes to G^^^ {x, x) and it is clear that G^^^ (a) and its derivatives 
diverge as cr„ tends to zero. These paths fall into classes Cn, which have n points on the 
boundaries. The length of the geodesic path that starts from x and arrives at x after n 
reflections from the boundary is shown by /„(x, x). These paths are closed in the limit 
X —>■ X. By imposing the above boundary conditions on the surface located at x^ = 
normal to the closed direction, the geodesic paths will be reflected on arriving at the surface. 
Dirichlet boundary condition on the plates is implemented by the reflection factor (-1). The 
generic closed path starting at x, reaches x„ after n-reflections from the boundary, where in 
the Cartesian coordinates x„ is given by 



n 1 2 \ 

y-'^ni •^ni •^ni •^n) 



(x°, 3x^ + {rio — 1) a, x^, x^) 
(x°, —x^ + iuo + 1) a, x^, x^ 



(11) 



in which rig and are even and odd integers respectively. It is convenient to separate the 
geodesic paths into odd and even paths, according to their number of reflections from the 
plates. Some of these paths are shown in Fig. 1. As it is seen, the two paths 2 and 2' with 
two reflections are distinct, so the multiplicity factor, M„ , is two for even paths, and one 
for odd ones. It is shown that odd and even paths differ in their contribution to the Casimir 
effect The contribution of the even paths in the vacuum energy density (0 |Too| 0) being 
finite, results in the Casimir force, while the odd paths lead to a divergent contribution. 
The half squared geodesic distance for a path with Ue = 2r refiections is 



X 



2raY - (x^ - x^)^ - (x^ 



X 



)• 



(12) 




X =a 



X' = 



r 



2' 



3' 



4' 



FIG. 1: Different geodesic closed paths. The initial and final points, which coincide, have been 
separated so that the paths could be identified. 

and its contribution to (0 |Too| 0) is 

1 



ne=2r 



limM 



2r 



(13) 



Since V^cr„^ = — V^(T„^ the contribution from even paths is given by 

oo 

{^\Too\0),,en = -lim VVoVoG«(a2,) 



m 



r=l 

4 oo 

hm > 



27r2 X 



K2 (A 



2r 2 

— m 



3- lVo(T2rJ 



(14) 



with = ml2r = 2mar. Note that the paths in class Co (r = 0) are excluded in the 
above summation as they lead to divergences. Since the interval between two endpoints is 
spacelike, that is cr„ > 0, only the first term in (6) concerned us in arriving at (fT^ . In 
coincidence limit, x ^ x, we have 



(0|Too|0)^ 



671 



2 °o 
2 



m" \ ^ K2 {2mra) 



r=l 



j,2g^2 



(15) 



Energy per unit area is obtained by integrating over the distance between the plates and 



is given by 



K2{2mra) 



This is a rapidly convergent series where the main 



contribution comes from the first term. Note that this result is valid for an arbitrary value 
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of In the limit ma <S 1, we obtain 

1 ^1 

-7r2^3 ~4 



r=l 



C(4) 



TV 



1440a3' 



(16) 



In the hmit m — > 0, the first term includes QO/tt^ fa 0.92 of the whole series. Next consider 
the odd paths. There are two famihes in each 6*2^+1 class illustrated in Fig. 1. One family 
begins its path with the first reflection from the lower plate while the primed family starts 
with reflections from the upper one. Half squared geodesic length for a path in class C2r+i 
is given by 



0"2r+l 



2r+l 



1 

2 



{x^ - 3x^ - 2raf 
{x^ + x^ -2{r + l) af 



( 



- x'^Y - (x^ - x^Y 



(17) 



For the tangent vectors we have 



In the coincidence limit, the lengths of the paths in class C2r+i range from 2ra to 2 (r + 1) a 
as varies from zero to a. The paths with Uo — 2r + 1 reflections contribute in energy 
density as 

(0|Too|0)2,+i 

= ^ lim {[-VoVo - (4^ - 1) ViVi] G'(i) ((72.+1) + [-V0V0 + (4^ - 1) ViVi] G^'^ } 



m 



47r^ x—yx 



I zr-|-i 

_^2^^3i_Wi) [(Voa2.+i)^ - (4, - 1) (Via2.+i)^] 
-m'^^^P^ Vyo^^.f - (4. - 1) 



(19) 



where ^2r+i = i^hr+i and ^2^+1 = i^l'2r+i- coincidence limit, we have 

, , , , r, , A'2 (2/// \.v^ + /-ol) A'2 (2//; |,r^- - (/• + 1) a 

Too 2, , , = — <^ 1 - 20 ^ ^ + 2<?^^ ^ ^ ^ 

\ I 001 /2.+1 27r2\ ^' 2w?{x^ + raf 2w? {x^ - {r + I) af 

r i^3(2m|x^ +ra|) _ K^{2m\x^ - (r + l)a|) " " 
^ ^ [ 2m (xi + ra) 2m {x^ - (r + 1) a) 



(20) 
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In the above expression we have set z = + ra\ or — (r + 1) a|. So the contribution of 
eq. to Casimir energy per unit area is 

(r+l)a 

/ {0\Too\0)dz. (21) 



Using the behavior of the function Ki, (z) for small values of its argument and the positive 
real part of u [l^ 

\im K,{z)^ It {u)(^Y\ (22) 



it is easily seen that the only divergent contribution comes from the paths in class Ci (r = 0) 
at the lower limit of ^2\\ji. It does not depend on relative variations of boundary surfaces 
and does not contribute to any forces. As discussed in the introduction this introduces 
a cutoff at a distance scale of A^^, where A is a frequency cutoff. In the regularization 
scheme used, this means that the minimum geodesic distance is chosen iohe x — x = hr^. 
Inserting A~^ as the lower limit for r = and summing over r, we obtain the contribution 

oo 

of odd paths eodd = f {0\Toq\0) dz. We note that the integrand leads to a quartic, A^, 
divergence which results in a cubic, A^, divergent Casimir energy. This contribution displays 
the expected cubic divergence for a scalar field obeying a Dirichlet boundary condition. Its 
finite contribution is canceled by the higher odd paths 

a 2a (r+l)a 

eodd = j {0 iTool 0)dz + J (0 \Too\0) dz + . . . + j (0 |Too| dz + . . . . (23) 

a ra 

Therefore the Casimir energy is given by the familiar result of eq. (fT^ . 



B. Permeable plates 

For an infinitely permeable boundary the use of closed paths renders the Casimir energy 
analysis transparent and rather trivial as we only need to replace the boundary conditions 
in (jH)) with, 

$ (x°,0,x2,x3) = 9^.1$ (x°,x\a;^x3) 1^1^^ = 0. (24) 

Note that the reflection factor will be (+1) instead of (-1). For the closed path of length 
= 2ra, the reflection factor (— 1)'^(+1)' = (—1)'^ must be inserted in the summation 
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(fl3j). so that the result is 



(Olrc^lO)^!^!: '-^''^/^""' . (25) 

which is similar to the case of twisted scalar fields discussed in Q|. In the massless limit, it 
leads to the energy per unit area 

r=l 

IV. CASIMIR ENERGY IN SPACETIMES WITH NON-TRIVIAL TOPOLOGY 

If we choose a periodic condition in direction with the periodicity equal to b, we get a 
cylindrical topology = R x S x R with the conditions 

$ 0, x^, x^) = $ a, x^, x^) = 

(J) /V»0 /vi2 /y»3\ (Jj / /v»0 /v^l /y»2 I qJ-^ -T*^^ 

and a non-zero vacuum energy is expected. The above introduced closed geodesic paths 
start at x and reach Xns such that 



(x° 



_i _2 ~S \ 
ns5 -^ns' ■^nsi •^ns) 



X"*^ + ngO, X^ + s6, X^) 

(x°,3xi + K-l)a,x2 + s6,x3) (28) 
(x°, — X"^ + ijio + 1) a, x^ + s6, x^) 



and n, s are integers. These paths can be classified into Cns classes, with n points on the 
plates and |s|-circulations in the x^ direction. Negative values of s refer to those paths 
that circle in the opposite direction. We see that G'(^) (x, x) is given by G'^'^'^ (x, x) = 
G*^^^ (x, Xns)- The tangent vector to the geodesic paths at x and Xns are related through the 

71, S 

equations V^o"„g = — V^cr„^ and (fTH|) . The lengths of the paths in the class C2r,s is fixed by 

r 2 2 ~i "~" r 2 2~i "~" 

the pair (r, s) while that of C2r+i,s varies from [(2ra) + {sh) J ^ to [(2 (r + 1) a) + {sh) J ^ . 

The zero path length belongs to the high energy fields in the class Co,o- Like the previous 

example, in the space integration we have 

e2r+M= j (0|Too|0)rf2, (29) 
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and for the class Ci_o, the distance scale enters at the lower limit and the finite contri- 
butions are canceled by the higher C2r+i,s classes. We get the finite results from the paths 
in C2r,s with (r, s) 7^ (0, 0). In this topology, the Casimir energy per unit area is 

m^a Ko (zr 



E (30) 



47r^ ^ z 

r,s=-oo,{r,s)j^{0,0) 



where [{2raf + {shfy. Our result in ^ would be in agreement with that ob- 

tained in for S = 5^ X i?, if the periodic length in direction is chosen to be 2a. Later we 
will use the analogy between the above Casimir energy calaculation in a nontrivial toplogy 
and that at finite temperature, to derive the free energy of a cavity at finite temperature. 
Using similar arguments for a spacetime with periodicities h and c in the and directions 
and the boundary conditions 

$ (x°,0,x^,x^) = $ (x°,a,x^,x^) = , , 

(31) 

$ x^, x^) = $ x\ + sb, x^ + Ic) , 

we obtain 

, _ K2 {Zrsl) 



r,s,Z=-oo,(r,s, 0^(0,0,0) 

where Zrsi = m [(2ra)^ + + (/c)^] ^. It can be seen that Zrs and Zrsi are nothing but 
the mass times the length of classical periodic path. A typical non-planar path starting in 
a given direction includes |2r|-refiections from the surfaces = 0, a, |s|-circulations in the 
x^ direction and |/|-circulations in the x^ direction. Again negative values of s and / refer to 
those paths that circle in the opposite direction. A factor of two in 2ra term is due to the 
Dirichlet boundary condition in the x^ direction. If the periodic condition in x^ direction is 
replaced with a Dirichlet boundary condition, we get the similar Casimir energy expression 
as in eq. ()32|) but now with Zrsi = rn [(2ra)^ + (2s6)^ + (/c)^] ^ . The Neuman boundary 
condition that the normal derivative of the field, rather than the field itself vanishes at the 
surface leads to the insertion of either (— 1)*^ or (—1)'^ for each permeable wall in x^ or x^ 
direction respectively. 

Since the fields at x and x„ {n a family index) are related through $ (x) = (—1)"$ (x„), 
applying the optical approach to twisted fields, different paths appear in eq. (fTn|l with non- 
equal weights. So a twist in any specific direction inserts a factor of (— 1)*^, where n is the 
number of circulations in that direction. 
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A. Non-orientable manifolds 



The manifolds we have examined so far are orientable. Here we examine the effect of 
spacetime manifold orientation on the closed geodesic paths and consequently on the vac- 
uum polarization. Since frames could only be defined locally and transport around closed 
spatial directions may change the handedness of local frames, it is expected that the non- 
orientability of a manifold could affect the Casimir energy calculations. A four dimensional 
Mobius strip xT, = x M"^ x is a non-orientable spacetime manifold which can be con- 
structed by the identification of the points (x", x^, x^, x^) and = {x^, x^ + na, (—1)" x^, x^) 
in Minkowski spacetime. So for a scalar field in this manifold we have 



$ x^ x^) = $ (x°, x^ + na, (-1)" x^ x^) . 



(33) 



Different paths in class C„, circulate |n|-times around x^ direction and reach x„ with the 
half squared geodesic lengths equal to, 



0"2r 



1 

2 



x°-x°)' 



(x^ — x^ — 2ra)^ + (x^ — x^)^ i , „ ^„,2 
(xi - xi - (2r + 1) a)' + (x2 + x2)' ' 



(34) 



The lengths of the paths in classes and C2r+i are given by \2ra\ and 

(2r + 1)^ + (2x^)^ ^ respectively. The zero-length paths fall into Cq class and so the 
finite contribution to the Casimir energy is due to the other classes. Now using the facts 
that, 

V^0"2r = -V^Cr2r 'V^a2r+1 = -V^(T2r+l + 25^^2 Via2r+1 , (35) 

we have 

(0 iTool 0) = -- lim V [2VoVoG'(i) (^2.) + [VqVo - (4^ - 1) V2V2] (^2^+1 

r=l 

_ nv' K2 {A2r) I 



27r' 



2 (1 - 2^) 



K2{A 



2r+lj 



(A 



(4(^-1) (2x2m)' 



2r-+lj 



2r-+lj 



2r+lj 



(36) 



where A2r = 2mra and A2r+i = rn (2r + 1)^ + (2x^)^ ^. This expression has been 
obtained in 13] using the method of images. 
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B. Mass and free energy expansion 



The above classical closed paths formalism could also be employed to calculate the free 
energy of the cavity discussed in the previous section at finite temperature jl^. To do so 
one could use the well known fact in the field theory at finite tempertaure that a system 
in three dimensional space in equilibrium at temperature T could be treated as a system 
in a four- dimensional Euclidean space whose fourth dimension is compactified to a circle of 
circumference It = This amounts to replacing in (jHUj) the contribution of each closed 
path of class C„ at T = by the sum of contributions of closed paths, now with added 
|s|-circulations in the extra fourth dimension. This extends the path lengths from In to 
C = [^n + ("5^r)^] ^ • So, using (jHn|l the free energy per unit area F is given by 

n=l s=-oo V " / 

Note the difference between eqs. (jHUj) and (jHTf) . When the boundary conditions ^Tl\ are 
applied, the eigenmodes will be labeled with two quantum numbers and consequently closed 
paths will be classified in classes labeled with two indices and, as before, zero length paths 
should be ignored. On the other hand for T > 0, the eigenfrequencies of the cavity change 
from ujn to uj'^ while the quantum numbers of a mode are conserved J^J . So only the paths 
in classes C^^o "will contribute while their lengths are given by instead of /„. In other 
words, there are no contributions in (jHTjl due to classes C„=o,s^o and the finite temperature 
does not induce any new ultraviolet divergences. Using (j37|l the free energy density can be 
written in the following form 



^4 ^ ^ (/T) 



71=1 S = -00 V 71 ) 

where its high and low temperature expansions can be obtained by introducing the variables 
T = 2akBT and a = m^j2a/kBT. may be rewritten as = a (^n^r + ^, and the 

summation ^ j, a is found by using the identity 



s=— 00 



/ (n) = —the sum of residues of vr cot (nz) f (z) at all poles of / (2;) . 

71= — 00 

In our case there are two poles at s = =i=mr and the residue at s = inr is given by 

TT^ r ,9, vrcothfTr™) vr coth (tt™) , , 

[1 - coth^ (vrrn)] - + ' > . (39) 
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As this is symmetric under n 



-n, the residue at s 



-iriT is the same and we obtain 



K2 m 



TT 



[1 - coth^ (vrrn)] + ^i"^- (40) 



2a^n 



It could be seen that after multiphcatiom by m"^/27r^ and summation over n, only the last 
term depends on the mass field and its contribution reduces the energy density. It is also 
seen that the mass term is completely separable and divergent in large r limit. It also 

oo 

behaves like ^ which introduces a new divergence. For massless fields and large r limit 

n=l 

it corresponds to 



ill 



yx2 



•2nnTk 



+ 



k=l 



TT 



1 + 2E' 



-2iTTnk 



k=l 



(41) 



and the free energy density is given by 



27r2 



TT 



27T 



aS) + — J2{l + 27rTk)e~ 



■2'KTk 



k=l 



J2 n 

n/ k=integer>0 



87ra3 



(42) 



where ( (3) = n ^ and a = a/m = ^J^ajk^T ■ Free energy at high temperature (fixed 

n=l 

a) converges like ^ \. Although its convergence is slower than that in the case at T = 0, 

n=l " 

the contribution of paths with small number of refiections is still a good apprroximation. 
To study the small r (low temperature) behaviour of the free energy density we first note 
a refiection symmetry r 1/rinjF — ^ ^^^""2^'^ which relates the high and low 

n=l 

temperature regimes. So with r replaced by 1/r in this expression, it is suitably adapted for 
obtaining the small r limit. Doing so we get another expression for eq. (|38|) in the following 
form 



m 
2^ 



E 

n=l 



[l - coth^ (vrn/r)] 



TTT COth {-KTi/t) 



7rcoth(7™/r) ^ K2 {an^) K2 {an/y/r) 



The last three terms show the mass dependence. Although from the previous studies [3]- 



it is known that the contribution of massive fields can be neglected relative to their 
massless counterpart, here infinities show up in low as well as in high temperature regimes 
whenever the fields are massive. The zero-point energy of a confined massive scalar field 
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in two dimensions is studied in [15| (see also Q]) where the mass dependent logarithmic 
divergence is absorbed into a redefinition of the so called bag constant. But in our case it 
is a problematic divergence as it can not be absorbed into a constant. In the massless case 
and for r ^ 1 we have 



TT 



n^ksT f , as) {ksT f , {ksT) 



+ 0{e 



1440a4 



-2-K/akBT 



90 



) =^0 



+ 



720 / T 



TT" 



1 + 27r- 



T 

Teff 
eff 



Ana 

4 

+ 



+ 



na 



TT 



aksT 



-n/aksT 



45C (3) 



eff 



T 



An 



(44) 

Here JFg is the zero-temperature Casimir energy density and Te// = is the effective 
temperature. The results (1421) a nd (1441) c an be compared with the high and low tempera- 
ture expansions obtained in [l^ and |l2|. Presence of the term proportional to in low 
T-expansion is an expected result when we note that it has the same form as the free en- 
ergy density of black body radiation of the same temperature in the same volume. The 
appeareance of this term originates from the fact that, as in the case of zero temperature, in 
renormalization of the free energy of the cavity we need to subtract from it the free energy 
of the unbounded empty space at the given temperature in the same volume (For a detailed 
discussion on this refer to and references therein). 



V. OPTICAL APPROACH AND SPINOR FIELDS 



In j8| we have discussed different spin structures, twisted and untwisted, on spacetime 
manifolds with non-trivial topologies by choosing different local frames. Analogous to eq. 
()1()|) . for spinor fields we have 

(0 IT'^I 0) = 4 lim Tr [(7^^^^ + 7"^ V'') - (7'^V'^ + Y^^] 7"V^ V G« (x, x„). (45) 

n 

For untwisted spinor fields and orientable manifolds, eq. (jl^ reduces to (0|T^i,|0) = 
lim;j^a; V^V^G'-^'' (x, x). Making use of the optical approach, different closed paths ap- 
pear with the same weight, except the zero length ones which should be discarded. Like 
the scalars, for twisted spinors the factor of (—1)" enters in eq. (j45|) to show the twist 
of the frame in the direction with n-circulation. In a non-orientable manifold such as the 
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Mobius strip, we will see that limited number of paths contribute due to the fact that the 
local frames at x and x„ (which are related through a Lorentz transformation) are not the 
same. The spinor \1/ (x„) is related to \E' (x) through \1/ (x„) = S'\l/ (x), where 5 is a suitable 
representation of SL (2, C). This translates into the Casimir energy expansion in the form 
of different weights for different closed paths. In Mobius strip, the local frames at x and x„ 
are related through parity and S = (7°7^7^)" such that S"^ = I and therefore, 

oo 

Gilia) = Yl KV2n) + G«(a2n+i)(7°7'7^)] • (46) 

n=—oo,nj^O 

Since the trace of an odd numbers of 7-matrices vanishes, taking the trace, reveals that the 
second sum in eq. fl46p makes no contribution and we immediately gain 

(oK„io) = -^E%^. m 

^ ~l {2mra) 

which, as expected, is minus the twice of the first term in eq. (jHBj) . 

VI. CONCLUSIONS 

Following Casimir's original prediction for the force between grounded conducting plates, 
many other calculations have been carried out for different variations in the conductor geom- 



etry or background 



, 3| and many exact and approximate techniques have been proposed 
to calculate this energy. Here as the main objective, employing the recently proposed optical 
approach, we calculated the Casimir energy for massive fields in a cylindrical topology as 
well as on a Mobius strip as a non-orientable manifold. The simplification provided by this 
approach allows the extension of the domain of applicability to the calculation of the reg- 
ularized vacuum energy of a quantum field subjected to various static external conditions. 
Applying this approach to systems at finite temperature, the finite temperature effect on 
the free energy of a cavity is compared with that of a periodic boundary condition at zero 
temperature. The same approach has also been applied to calculate the Casimir energy for 
spinors and twisted fields, where it is shown that different paths contribute with non-equal 
weight factors. Using the same approach we have recovered the well known fact that the 
dependence of free energy on mass at high and low temperatures may be regarded as prob- 
lematic since the key feature here is the presence of a divergent term which is absent in the 
massless case. As a consequence the contribution of massive fields in these regimes not only 
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does not tend to damp relative to their massless counterparts (as in the T = case) but 
also diverges. 
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